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ABSTRACT 



We present a systematic and physically motivated char- 
acterization of incoherent or coherent electric and mag- 
netic fields, as measured for instance by the low fre- 
quency receiver on-board the Solar Orbiter spacecraft. 
The characterization utilizes the 36 auto/cross correla- 
tions of the 3 + 3 complex Cartesian components of the 
electric and magnetic fields; hence, they are second or- 
der in the field strengths and so have physical dimension 
energy density. Although such 6x6 correlation matri- 
ces have been successfully employed on previous space 
missions, they are not physical quantities; because they 
are not manifestly space-time tensors. In this paper we 
propose a systematic representation of the 36 degrees-of- 
freedom of partially coherent electromagnetic fields as a 
set of manifestly covariant space-time tensors, which we 
call the Canonical Electromagnetic Observables (CEO). 
As an example, we apply this formalism to analyze real 
data from a chorus emission in the mid-latitude magne- 
tosphere, as registered by the STAFF-SA instrument on 
board the Cluster-II spacecraft. We find that the CEO 
analysis increases the amount of information that can be 
extracted from the STAFF-SA dataset; for instance, the 
reactive energy flux density, which is one of the CEO pa- 
rameters, identifies the source region of electromagnetic 
emissions more directly than the active energy (Poynting) 
flux density alone. 

Key words: electromagnetic, observables, correlation, 
space-time, covariant, irreducible, tensor, Cluster 



1. BACKGROUND 



When analyzing time varying electric and magnetic vec- 
tor field data from spacecraft, it is common to construct a 
6x6 matrix from the complex sixtor (E, B). In a Carte- 



sian coordinate system, this matrix can be written 
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This electromagnetic (EM) sixtor matrix has in various 
guises, such as Wave-distribution functions (WDF) ([l|) 
and so on, been useful in the analysis of EM vector field 
data from spacecraft, such as on the Cluster and Polar 
missions. Although this matrix-form description of the 
second order properties of the EM fields can in some in- 
stances be convenient and intuitive, it unfortunately has 
no real motivation in physics, since physical quantities 
such as EM fields are ultimately not sixtors. 

The EM sixtor matrix lists all possible auto and cross cor- 
relations of E(x, t) and B(x, t) and, hence, contains the 
complete information of the second order properties of 
EM waves. It is clear that one obtains the EM energy 
densitjQ |Ep + |Bp by taking the trace of the EM six- 
tor matrix but it is not clear what other EM quantities 
can be extracted and how this extraction should be per- 
formed in general. Some EM quantities are obvious and 
can be picked out by hand, such as the Lagrangian den- 
sity |Ep - I Bp, the Poynting flux density K{E x B}, 
etc; B denotes the complex conjugate of B, etc. Other 
EM quantities are not so easily identified. Furthermore, 
as it stands here, the EM sixtor matrix is a mixture of 
both scalar and pseudo scalar components. This is due 
to E being a proper (polar) vector and B being a pseudo 
(axial) vector. This could be remedied by for instance us- 
ing the sixtor (E, ±iB) but then again, it is unclear what 
sign to use for the imaginary unit, i. Regardless of how 
the sign of i is chosen one would need to redefine the EM 
quantities. The trace of the modified proper EM sixtor 
matrix would in this case correspond to the Lagrangian 
density rather than the EM energy density, which in turn 
would have to be redefined. Other EM quantities would 
also need to be redefined in non-standard ways. 

' Throughout the paper, when second order field quantities are dis- 
cussed we have chosen a normalization such that the speed of light is 
set to unity. 



Perhaps even more important for space borne observa- 
tions: the EM sixtor matrix is not covariant according 
to the requirements of special relativity. Spacecraft are 
constantly moving and often spinning observation plat- 
forms. EM wave measurements become Doppler shifted 
and data must often be_"despun". For the 3x3 sub- 
matrices, E ® E, E ® B, B (g) E, and B (g) B, where (g) 
denote the direct product, despinning is straight forward 
by applying rotation matrices R from left and from 
right, e.q. E' (g) E' = RE (g) ER'^. To rotate the full EM 
sixtor matrix, similar operations must be performed four 
times. This is awkward and the resulting 6x6 matrix is 
still not covariant. For EM wave measurements in space 
plasma the last remark can be crucial. 

A Lorentz boost is the translation from one Lorentz frame 
to another one moving at velocity v. A Lorentz boost 
does not necessarily imply relativistic speeds, which is 
a common misconception; and therefore it do not by it- 
self preclude what is typically associated with relativis- 
tic effects. It is simply a quite general recipe to make 
two different observers agree on a physical observation. 
The Lorentz boost of the EM field vectors can be written 
E' = 7(E+vxB)andB' 7(B-vxE/c2), where cis 
the speed of light and 7=1/ -^/l — jc? . As a matter of 
fact, the Lorentz boost is the essence of the well-known 
frozen-in field line theorerrQ from magnetohydrodynam- 
ics (MHD); a theory which is commonly used to model 
the solar wind plasma. In a plasma, relativity comes into 
play at very a fundamental level since the electromagnetic 
(Lorentz) force dominates the vast majority of all plasma 
interactions. 

Another example illustrates the problem to separate time 
(frequency) and space (wave vector) in EM wave obser- 
vations on board a spacecraft. Assume that we observe a 
wave mode which is described by an angular frequency 
w and wave vector k. We can write this as a 4-vector 
(w, ck). Let's make a Lorentz boost in the v direction: 



ijj' = 7(0; — k • v) 
[7-1 



ck' = ck 



V 



(ck ■ v) — 70; 



(1) 
(2) 



What happens now for a stationary (DC) field structure 
moving with the solar wind plasma? We then have cj = 
and |ck| 7^ 0. For a satellite moving with velocity v 
relative to the DC field structure, it is justified to set 7 « 1 
(the solar wind speed seldom reaches more than 900 km/s 
and using this value we obtain 7 w 1.0000045 > 1); Eqs. 
([T]i and (|2|i are then reduced to 



ck' 



-k - V 

ck 



(3) 
(4) 



We can see that the DC field structure is not Lorentz con- 
tracted appreciably at this low velocity, k' w k. How- 
ever, there is a dramatic change in the observed fre- 
quency, which for a head-on encounter with the structure 



^If E + V X B = in a plasma, the magnetic field lines change as 
though they are convected with velocity v, i.e., they are frozen to the 
plasma flow. This is the frozen-in field hne theorem of ideal MHD. 



is registered as uj' k, kv rather than zero. The observed 
frequency is proportional to the dimension of the struc- 
ture, which we take to be in the order of one wavelength, 
A = 27r/fc. Taking v = 900 km/s a 900 km DC field 
structure would now register as 1 Hz, a 90 km structure 
as 10 Hz, and a 9 km structure as 100 Hz, etc. 

These simple examples clearly show that a space-time 
(covariant) description is necessary even if 7 > 1. The 
frequency (time) and the associated wave vector (space) 
can not be treated separately but must be considered to- 
gether, as a space-time 4-tensor 

The Maxwell equations are inherently relativistic and can 
easily be put into a covariant form using 4-tensors. From 
a theoretical point of view, this fact alone provides a very 
good argument why one should try to express also the 
second order properties of the EM fields using a covari- 
ant formalism. This was recently carried out by the au- 
thors and published in a recent paper In this pa- 
per we introduced a complete set of space-time tensors, 
which can fully describe the second order properties of 
EM waves. We call this set of tensors the Canonical Elec- 
tromagnetic Observables (CEO); in analogy with Wolf's 
analysis of the Stokes parameters We suggest that 
the CEO could be used as an alternative to the EM six- 
tor matrix. Not only are the CEO covariant, but they 
are all real valued and provide a useful decomposition of 
the sixtor matrix into convenient physical quantities, es- 
pecially in the three-dimensional (3D), so-called scalar- 
vector-tensor (SVT) classification; see section l272l The 
CEO have all dimension energy density but have various 
physical interpretations as will be discussed in what fol- 
lows. 



2. CANONICAL ELECTROMAGNETIC OB- 
SERVABLES 



The CEO set was derived from the complex Maxwell 
field strength F"^. Other possibilities, such as using the 
4-potential or using a spinor formalism (4), were con- 
sidered but discarded due to their lack of physical con- 
tent. The 4-potential is not directly measurable and it is 
furthermore gauge dependent. Spinor formalism has been 
proved possible to use Q but we believe the space-time 
tensor formalism to be more intuitive and convenient to 
use. 

In the quantum theory of light, observables of an EM field 
are ultimately constructed from a complex field strength; 
see (111). The simplest of these observables are sesquilinar- 
quadratic (Hermitian quadratic) in F^^, i.e., they are 
functions of the components of F<^0pi^^ which is a 4- 
tensor of rank four We showed that it was possible to de- 
compose F""^ F^^ into a unique set of tensors, the CEO, 
which are real irreducible under the full Lorentz group. 
We shall not repeat the derivation here but will instead 
discuss the space-time (4-tensor) and three-dimensional 
(3-tensor) representations of the CEO. 
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Table 1. CEO in space-time classification, i.e., 4-tensor 
notation: 1 + 1 + 9 + 9 + 6 + 10 = 36 observables. 



2.1. Fundamental space-time representation 

In terms of the Maxwell field strength F'^^, the CEO 
are organized in the six real irreducible 4-tensors 
C+, C_, Q"^, T"'^, C/"^, and W"'^^^ . This is the funda- 
mental space-time representation of the CEO; their prop- 
erties are listed in Table l2.l] 

The CEO 4-tensors are defined as follows: the two scalars 
are the vacuum proper- and pseudo-Lagrangians, 

C+ := {F.^pF'^f' ~ *F^p*F^'') /2, (5) 
C_ := {F^p*F"^ + *F^pF^P) /2, (6) 

respectively, where we have used the dual of F°'^ defined 
as 

*p^p _ ie"'5''^F^5 = ^e'^^^sF'"'- (7) 

The three second rank tensors consist of the two symmet- 
ric tensors 

._ (^po^^pf^P ^ *pypt,l3'^ /2, (8) 
jjat3 _j (^pypt.13 _ *j7'"^^M/3) /2, (9) 

and the antisymmetric tensor 

ga/3 _j (^po^^pt^a _ *po.^*pf^P _ /2. (10) 

The symmetric second rank tensor T"'^ is the well-known 
EM energy-stress tensor, which contains the total energy, 
flux (Poynting vector), and stress (Maxwell stress tensor) 
densities. The other two second rank tensors, U"'^ and 
respectively, are less well-known. The symmetric 
1/°"^ tensor is similar to T"'^ in that it contains active 
energy densities but in U"^ these densities are weighted 
and depend on the the handedness (spin, helicity, polar- 
ization, chirality) of the EM field. Therefore, we have 
chosen to call them "handed" energy densities. The anti- 
symmetric tensor Q"^ on the other hand is very different 
in that it only contains reactive energy densities, which 
are both total (imaginary part of the complex Poynting 
vector) and handed. 

The fourth rank tensor is 

- '^C+jf^^rf^f^ - -C-e°'f^^^ (11) 
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Table 2. CEO in scalar-vector-tensor (SVT) classifica- 
tion, i.e., 3-tensor notation: 4x(l + 3 + 5) = 36 observ- 
ables. 



where the square brackets denotes antisymmetriza- 
tion over the enclosed indices, e.g., T"-^^g"'^^ = 
1 (jaSg^p _ T^-tgSP^^ and nested brackets are not op- 
erated on by enclosing brackets, e.g., T^°'^^ g'l^l^^ = 

1 ^paSg-yP „ patgSp _ pPSgfa _^ pPfgSa^ ^ jj f^jfjUg 

the symmetries W"^'^^ = W^""'^ = W"^^'^ = W^^"^ 
and VV^^f^^^^ = 0. 

This real irreducible rank four tensor, Eq. (fTTT l. was dis- 
covered by u^ and published in (2), and is still under 
investigation; it is an extremely interesting geometrical 
object, having a structure identical to the Weyl tensor in 
general relativity; see (@). We have found that it contains 
a four-dimensional generalization of the Stokes param- 
eters, as will be demonstrated in section [5] for the two- 
dimensional (2D) case. It contains both reactive total and 
reactive handed energy densities. 



2.2. Three-dimensional representation 

The fundamental space-time 4-tensor CEO can be written 
in terms of the three-dimensional E and B vectors, i.e., 

3- tensors. This is convenient because it allows us to use 
intuitive physical quantities. To systematize the 3D rep- 
resentation of the CEO, we will use a physical classifica- 
tion where we organize the CEO into four groups, which 
have been introduced briefly in the previous section: the 
(active) total, (active) handed, reactive total, and reactive 
handed CEO parameter groups, respectively. In addition, 
we will use a coordinate-free 3D formalism and classify 
the CEO parameters according to rank, i.e., as scalars, 3- 
vectors, and rank two 3-tensors (SVT classification). The 
3D CEO are listed in Table IZ2] The CEO 3-tensors are 
defined as follows. 

The "total" parameters are: 

u =r°" = (|E|2 + |B|2) /2 (12) 
N =r'" = 5R{E X B} (13) 
M =r'^ = ulg - SR{E®E + B®B} (14) 

where I3 is the identity matrix in three dimensions. This 
is the 3D representation of the well-known energy-stress 

4- tensor T°''^, defined by Eq. dD. 



'To the best of our knowledge, the W"^^'^ tensor has never before 
been pubhshed in the hterature. 



The "handed" parameters are: 
h^U°° = 3{E-B} 



(15) 



S =[/* = ---3{(E X E + B X B)} (16) 



C =^U'^ = /il3 - 3 {E (g) B - B ® E} 



(17) 



This is the 3D representation of the handed energy-stress 
4-tensor J7"'^, defined by Eq. (|9]l. 



The "reactive total" parameters are: 
Z =C+ = (|Ep - |Bp) 12 



{E X B} 



R =Q'0 

1 / „ _ 2 

X =W'°^° = - 5R{E®E-B®B}- -/I3 

2 V 3 



(18) 
(19) 

) 

(20) 



Contrary to the active, total and handed, parameter 
groups above, the reactive total parameter group have no 
single corresponding 4-tensor Instead it is composed of 
parts from three different CEO space-time tensors: the 
vacuum proper-Lagrangian defined by Eq. (|5]i, the re- 
active energy flux density from Eq. ( fTOb . and the gen- 
eralized Stokes parameters corresponding to the auto- 
correlated E and B fields from Eq. (fTTI) . 



The "reactive handed" parameters are: 
a=C- = -SR{E-B} 



(21) 



0=iei,Q'' = -^^^{(ExE-BxB)} (22) 



Y =l-eiiW'°''^ = i ( K {E B + B ® E} - ^alg 



(23) 

Also for this parameter group, there is no single corre- 
sponding 4-tensor. The reactive handed group is com- 
posed of parts from three CEO space-time tensors: the 
vacuum pseudo-Lagrangian, defined by Eq. the re- 
active handed energy flux density from Eq. ( fTOl ). and 
the generalized Stokes parameters corresponding to the 
cross-correlated E and B fields from Eq. (fTTI) . 



(2D) fields in this plane be denoted E2D := {Ex , Ey) and 
B2D := By), and define the scalar product between 
2D vectors as 



E2D • — ExBx + EyBy 

and the cross product as 

E2D X B2D = Ex By — EyBx 

and the direct product as 

T71 TD / ExBx Ex By 

E2D <&) rJ2D = p R PR 



(24) 



(25) 



(26) 



We will however not need to consider all the components 
of the 2D direct product since the 2-tensors we will con- 
sider are all symmetric and traceless. Hence, we only 
want the parameters which correspond to Pauli spin ma- 
trix components 



1 

1 

1 

-1 



(27) 
(28) 



The Pauli components can be extracted from a 2D matrix 
by matrix multiplying by a Pauli spin matrix and then 
taking the trace, that is 

Tr {(E2D <S> B2d) <Jx} = (E2D <E) B2d) • -Tx 

= ExBy + EyBx (29) 

where we have introduced the double scalar product, ax ■ 
■ay, see (0). 

We can derive a set of two-dimensional canonical elec- 
tromagnetic parameters from the full CEO by formally 
taking 



= B, 







(30) 



and discarding all the parameters that are identically zero. 
In this way we obtain the following set, which we write 
in the coordinate-free 2D formalism introduced above. 



3. CEO IN TWO DIMENSIONS 



Up until now we have assumed that all three Cartesian 
components of both the electric field, E, and the mag- 
netic field, B, are measured. One may ask what happens 
if some components are not measured; can all the 36 pa- 
rameters of the CEO be retained? Of course this is not 
possible, some information is certainly lost in this case, 
but what one can do is to construct a set of parameters 
analogous to CEO in two-dimensions. 

Assume that we can measure the electric field and the 
magnetic field in a plane which we can say is the xy- 
plane without loss of generality. Let the two-dimensional 



The "total" 2D parameters are: 



U2D 



E. 



2D I 



|B2d| 

2 , 



/2 



\Bx 



By\'\/2 (31) 



\Ex\ +\Ey\ -r 
iV3 = K{E2DXB2D} 

= 5? {ExBy — EyBx} (32) 
M^^ = SR{E2D ® E2D +B2D ® B2d} • 

= {\Ex\'' - \Ey\'' + \Bx\'' - \By\''^ 12 (33) 

= K {E2D ® E2D + B2D ® §20} • •cra,/2 

= ^{ExEy+BxBy} (34) 



The "handed" 2D parameters are: 

h2D = 3{E2D -BaD} 

= {E.,B, + EyBy] (35) 

5*^ = 3 {E2D X E2D + B2D X B2d} 

^^{E,Ey+B.,By] (36) 

= 3 {E2D ® B2D - B2D ® E2d} ■ -(7^2 

= ^{EA-EyBy] (37) 

= 3 {E2D ® B2D - B2D «) E2d} ' ■<Jx/'2. 

^Q{E,By+EyB,} (38) 

The "reactive total" 2D parameters are: 

l2D= (|E2d|'-|B2d|')/2 

= [\Exf + \Eyf-\Bxf~\Byf)/2 (39) 

= 3{E2D X B2d} 

^^{E,By-EyBx} (40) 

X^^ = {E2D ® E2D - B2D (» B2d} • •crj2 

= {\Ex\'-\Ey\'-\Bxf + \Byf)/2 (41) 

= 5R{E2D ®E2D -B2D ®B2d} ' -(7^/2 

= 5R{£;,4-B,Sy} (42) 
The "reactive handed" 2D parameters are: 

0,20 = 3? {E2D ■ B2d} 

= ^{EA + EyBy} (43) 

= 3 {E2D X E2D - B2D X B2d} 

= ^{E,Ey^ B,By} (44) 

Y^^ = ^ {E2D ® B2D + B2D ® E2d} • ■cr./2 
= 3fi{S,B,-SyBj,} (45) 

= ^ {E2D «) B2D + B2D ® E2d} • -(Tx/^ 

= ^{ExBy+EyBx} (46) 

We can associate names with these parameters as listed 
in Table [3] The first four parameters, which we call the 
"total" 2D CEO parameters are all well known. These 
parameters are also known by different names, e.g., the 
total energy flux is also known as the Poynting vector 
(z-component), and the total energy stress is known as 
the Maxwell stress tensor (difference of diagonal com- 
ponents and off-diagonal component). The remaining 
three sets of 2D CEO parameters are less well known. 
We will not be able to provide a full physical interpre- 
tation of each of these parameters; indeed their role in 
space plasma physics is yet to be fully explored. We will 
only mention that the "handed" parameters involve spin 
(helicity, chirality, polarization) weighted energy, i.e., 
the energy of the right-hand wave modes are weighted 
positively and the energy of left-hand wave modes are 
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0, 
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EM Stokes parameter Q cross-type 




EM Stokes parameter U cross-type 



Table 3. Naming scheme for the 2D electromagnetic 
second-order parameters. All parameters are implicitly 
densities and two-dimensional. 



weighted negatively, and these weighted energies are then 
added. Its flux corresponds to the concept of elliptic- 
ity and for the case of vacuum, it is numerically equiva- 
lent to Stokes V parameter. The reactive energy densities 
come in two groups: the "reactive total" and the "reac- 
tive handed" 2D CEO parameter groups. From the "re- 
active total" group, we now recognize the reactive energy 
flux density, as well as the EM Stokes Q and U parame- 
ters, which here are of the auto-type; the vacuum proper- 
Lagrangian needs no further introduction. The "reactive 
handed" group contain the handed counterparts of the re- 
active energy flux density and EM Stokes parameters, 
which here are of the cross-typer; the vacuum pseudo- 
Lagrangian is well-known. 



4. APPLICATION OF CEO TO CLUSTER DATA 



Let us demonstrate that the CEO parameters can easily 
be computed from actual data. Assuming that we have 
measurements from a vector magnetometer and an elec- 
tric field instrument, all that is required is to auto/cross- 
correlate all measured components and then form the ap- 
propriate linear combination introduced above. 

As an example we will consider the STAFF-SA dataset 
on the Cluster-II space-craft mission (H. The STAFF- 
SA instrument dlOl) is well suited for the CEO parame- 
ters since it outputs auto/cross-correlation of electric and 
magnetic field components; however as Cluster does not 
measure one of the electric field components (namely the 
component normal to the spin-plane of the space-craft) 
we can only use the 2D version of the CEO introduced in 
the previous section. 

For this particular example, we re-process the high-band 
part of STAFF-SA data from an event discussed in @ 
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Figure 1. Example dynamic spectra of the 16 normalized two-dimensional CEO parameters. The parameters were 
computed from STAFF-SA data from Cluster space-craft 2 using the ISDAT database system. The following normalization 
has been applied: each parameter has been divided by the total energy except the total energy itself. Thus all spectral 
values are in dimensionless unit except for the total energy. This Figure can be compared with Fig. 1 in d^. The 16 
parameters are subdivided into a) the total energy parameters, b) the handed energy parameters, c) the reactive total 
energy parameters, and d) the reactive handed energy parameters. Note that all the quantities are purely electromagnetic 
in origin and so do not refer to contributions from the plasma. 



from 2001-03-31 UT. In Fig. 1 of this paper. Parrot et 
al display certain parameters based on the STAFF-SA 
data computed using a numerical software package called 
PRASSADCO; see diH). The interesting feature of the 
2D CEO parameters is that they are the complete set of 
electromagnetic field observables in the spin-plane of the 
space-craft; and indeed, they use up all the parameters 
in the STAFF-SA dataset expect for the magnetic field in 
the spin direction. Each CEO is a distinct physical quan- 
tity and examination of the panels in Fig. [T]indicates that 
this is indeed the case, since besides showing a common 
chorus feature (the arch to the left in each panel) there are 
unique points in each of the panels. 

Besides being a complete description of the electro- 
magnetic observables, the fact that the CEO parameters 
are based on parameters that conform with the physics 
of space-time means that we can expect physical phe- 
nomenon to be measured properly. Seeing as how the 
CEO parameters have not been explicitly measured in the 
past, we can expect that their future use may lead to new 
physical insights, especially since several of the parame- 
ters are completely new to space-physics. As an example 
consider again the data shown in Fig [T] It is interesting 
to note that the reactive total energy flux is only signifi- 
cant close to the equator; this implies that the equator is 
the source region for the chorus events, since reactive en- 
ergy flux is typically large close to radiating objects due 
to large standing energy fields. One can also see a modu- 
lation at 2.5 kHz in the EM Stokes parameters. If this is 
a physical phenomenon it would be indicative of Faraday 
rotation. Also there seems to be frequency dispersion in 
the handed stress since its components changes sign with 
frequency. Finally, the handed energy clearly shows the 
handedness of the chorus emissions on its own, without 
recourse to the sign of the total energy flux. 



5. CONCLUSION 



The proposed CEO parameters conveniently organize the 
measurements of the full EM field. Furthermore, they are 
physically meaningful quantities, i.e. they 

• have conservation laws 

• transform as geometric (Minkowski space-time) ob- 
jects 

• are mathematically unique (they are irreducible ten- 
sors) 

• retain all information, i.e. nothing is lost (linear 
transformation back to full sixtor form exists) 

• enables considerable data reduction (through param- 
eter subset selection) 

• have clear despinning properties {e.g. scalar quanti- 
ties do not need despinning!) 

• are all real valued 



• provides useful decomposition of the 36 second or- 
der EM components into twelve 3 -tensor quantities 

• reveals some new physical parameters describing 
EM waves: opening for new physical insights. 
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